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. . . The gravitational self-force on a point particle moving in a vacuum background spacetime can be 

' expressed as an integral over the past worldline of the particle, the so-called tail term. In this paper, 

I we consider that piece of the self-force obtained by integrating over a portion of the past worldline 

■ that extends a proper time At into the past, provided that At does not extend beyond the normal 

(N . neighborhood of the particle. We express this "quasi-local" piece as a power series in the proper 

^ I time interval At. We argue from symmetries and dimensional considerations that the 0(At'') and 

^ , 0(At) terms in this power series must vanish, and compute the first two non-vanishing terms which 

' occur at O(At^) and O(At^). The coefficients in the expansion depend only on the particle's four 

velocity and on the Weyl tensor and its derivatives at the particle's location. The result may be 
useful as a foundation for a practical computational method for gravitational self-forces in the Kerr 
spacetime, in which the portion of the tail integral in the distant past is computed numerically from 
a mode sum decomposition. 



I. INTRODUCTION AND SUMMARY 



> 

o 
o 

One of the outstanding open problems of classical general relativity is the calculation of the gravitational self-force 
"nI . experienced by a massive particle moving in a curved background spacetime. Here by particle we do not mean a 
■ point particle, but rather an extended object whose internal structure has a negligible effect on its "center-of-mass" 
O . motion. Such a particle will not follow a geodesic of the background space-time, but rather a geodesic of the total 
'Y* spacetime, whose curvature reflects both the background and the mass/energy of the particle itself. However, if the 
mass of the particle is much less than the natural length-scale of the background (i.e. the square-root of the inverse 
of the curvature scale of the background), the deviation of the particle's trajectory from a background geodesic will 
J> , be small over time-scales that are less than the natural length-scale of the background. In this case, one can treat 
k>l ' difference between the total spacetime metric and the background spacetime metric to linear order as spin-2 field 

r> ' generated by the particle and living on the background spacetime. This field couples to the particle, and in this 
. picture, it is understood as causing a force, the self-force, which causes the particle to deviate from the background 
■ ■ ' geodesic. 

With this description, the gravitational self-force is analogous to the self-force experienced by a electrically charged 
point particle coupled to a Maxwell field in curved spacetime, or to the self-force experienced by a particle carrying 
a scalar charge coupled to a massless linear scalar field in curved spacetime. In each of these cases, the interaction 
of a particle with its own field alters its motion. In the flat spacetime limit, these forces reduce to the familiar 
radiation reaction forces, as they are only associated with accelerated motion and the resultant emission of radiation. 
In curved spacetimes, however, the notion of emission of radiation can, in general, be ambiguous. Furthermore, even in 
stationary spacetimes, the self-forces can have conservative terms, as well as dissipative pieces reminiscent of radiation 
reaction forces. Thus, in curved spacetimes, the notion of a radiation reaction force is replaced by the more general 
notion of a self-force. 

Recently, there has been renewed interest in finding the gravitational self-force experienced by a massive particle 
traveling in a curved background spacetime. The primary impetus for this interest is the imminent construction of 
space-based gravitational wave observatories, such as LISA 1]. One of the most interesting gravitational wave sources 
for these instruments will be compact solar mass objects inspiralling into black holes of 10"^ to 10^ solar masses out 
to Gpc distances^- The extraction of the maximal amount of information from such observations, however, can 
only be effected if accurate theoretical waveform templates exist. The calculation of such templates requires precise 
knowledge of the orbital evolution of the smaller object in the gravitational background of the massive black hole, 
which in turn requires an accurate calculation of the gravitational self-force. 

In flat spacetime, the radiation reaction force can be thought of as the recoil force experienced by the particle as 
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it emits radiation. Since the radiation carries momentum, and since the radiation field is not generally spherically 
symmetric, there is a non-zero net momentum transferred to the particle. This change of momentum corresponds to 
a force, which for an electrically charged particle in flat spacetime is given by the Abraham-Lorentz-Dirac formula 
1^1^. This explanation of the radiation reaction force is, of course, highly simplified. It is not straightforward to 
relate the force directly to the radiation — the force is local to the particle while radiation is a far-field phenomenon. 
Nonetheless, this simplified picture can provide intuitive insight into the phenomenon. 

In curved spacetimes, qualitatively new types of self interactions occur due to the failure of Huygens' principle (in 
its modern incarnation) to hold in most geometries ^. More pedantically, it is due to the fact that solutions to wave 
equations on most curved manifolds depend not only on Cauchy data directly intersected by the past light cone, but 
also on Cauchy data interior to that intersection. The portion of the field that propagates in the null directions along 
the characteristics is called the direct part. The portion of the field which propagates in the timelike directions in 
the interior of the light cone is called the tail. Clearly, the tail part of a particle's field can interact with the particle, 
leading to a contribution to the self-force. 

A general expression for the self-force on an electrically charged particle in a curved spacetime was obtained in the 
seminal paper by DeWitt and Brehme B ^. More recently, similar expressions have been obtained for the gravitational 
self-force by Mino, Sasaki and Tanaka |lOj | and by Quinn and Wald lllj, which we review in Sec. ITU below, and for 
the scalar self-force by Quinn 1 121. These results have resolved many of the issues of principle in computing self-forces 
in curved spacetime. See Ref. 'l3l for a detailed review of and simplified versions of these computations. 

However, for applications to gravitational wave observations, one needs to translate the formal expressions of Refs. 
[lOllllj into practical computational schemes for computing orbits of particles in the Kerr spacetime. The expressions 
for the self-force involve the retarded Green's function for the wave equation, and the standard method of computing 
this Green's function is to use a decomposition of the field into modes. This mode decomposition method combines 
together the tail and singular pieces of the fields in a manner that is difficult to disentangle, and it is the tail piece 
of the self-field that determines the self-force. Thus, the results of Refs. II in in do not directly give a simple 
method of computing self-forces in black hole spacetimes. 

Nonetheless, some progress has been made in calculating the self-force for particular particle trajectories in particular 
spacetimes. For specific geometries, and also in the weak-field approximation, it has been possible to compute the 
tail of the Green's function. Some classic results have been obtained for scalar or electric charges in static, radial, or 
circular trajectories about black holes or cosmic strings [lllllllllTlllllli.l20ll2ll2l.l23Ll23 |. In Schwarzschild, 
and for circular or equatorial orbits in Kerr, the time-averaged non-conservative (i.e. radiation reaction) contributions 
to the gravitational self-force may be deduced using ener gy a nd angular-momentum balance arguments involving the 
flux of radiation to infinity and down the black hole [2^ 12^ l27l l2a . |2^ Isol IsiL Is^ . Furthermore, in the weak-field, 
slow-motion limit of general relativity, one may use the post-Newtonian expansion to obtain the grav itational self-force 
^^Is^lssf. and the result agrees with that obtained by specializing the formal results of Refs. |ldll]| to weak fields 

While these results are encouraging, it is important and desirable to have a framework in which arbitrary motions 
in black hole spacetimes can be computed. Recently considerable progress has been made in developing practical 
computational schemes for obtaining self- forces Most of these schemes are based on computing the full retarded 
field, which is infinite on the particle's worldline, and regularizing it in some way to effect the subtraction of the 
direct part of the field, leaving the desired tail part. Barack and Ori have derived a mode sum regularization 
scheme |3^ IstI ,38.] that has been successfully applied in a number of cases 38, 39, 40, 41, 42, 43]. The regularization 
parameters for this scheme have been derived from the fundamental MSTQW equation of motion for Schwarzschild 
in Refs. 0,^^ and for Kerr in Ref. Mino, Nakano and Sasaki have developed two regularization schemes, one 
of which is a mode-by-mode regularization and the second of which they dub the "power expansion regularization" 
which involves a Post-Newtonian expansion of the Green's function 47, 48, "49]. Both their methods have been applied 
[set l5l|. Another scheme is that of Lousto, who used zeta- function regularization of modes for a radially infalling 



^ This statement is in fact a modification of a conjecture by Hadamard|^ that the only hyperboUc differential operators whose solutions 
obey Huygens' principle are conformally related to the the ordinary wave operator in Minkowski space with even numbers of spacetime 
dimensions. Hadamard originally formulated this conjecture with scalar operators in mind, but it has since been trivially extended 
to operators for fields with higher spin. Counter examples to Hadamard's conjecture have been found, originally in > 6 (but even) 
dimensions by StellmacherQ and later in the more physically relevant case of > 4 (even) dimensions by GiintherQ- Nonetheless, these 
counter examples are believed (but not yet proved) to be isolated cases. This modified Hadamard conjecture has been proved in broad 
classes of spacetimes (e.g. some algebraically special spacetimes) and work continues in this area. 

^ The expression obtained by DeWitt and Brehme is missing a term due to a trivial calculational error; sec Hobbs 9] for the correction. 

^ See the proceedings of the Capra Ranch meetings on radiation reaction at http : //www. Isc-group .phys .uwm. edu/' patrick/irelaiid99/, 
http : //www. tapir . caltech. edu/~ capraS/, http : //www. aei-potsdam.mpg.de/lousto/capra/, http : //cgwp.gravity .psu. edu/events/Capra5/capraE 
and http : //cgwa . phys . utb . edu/Event s/agendaView . phpTEvent ID3. 
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scalar particle in Schwarzschild 52]. 

A qualitatively different method for computing self- forces in black hole spacetimes has been suggested by Poisson 
and Wiseman 0|. It is based on a direct computation of the tail field, rather than a regularization of the total 
retarded field. The tail field can be expressed as an integral over the past worldline of the particle. The idea is to 
split this integral into two pieces, a piece that extends back into the past a proper time interval At, which we call the 
quasi-local piece, and the remainder of the integral. The second piece, from the distant past, can be computed using 
standard multipolar decomposition of the full retarded field; no difficulties involving divergences occur here, and thus 
no regularization is needed. The first, quasi- local piece can be computed approximately as a power series expansion 
in At. 

In this paper, we compute the expansion in powers of At of the quasi-local piece of the gravitational self force for 
an arbitrary vacuum spacetime, to the first two non-trivial orders in At. Our result is given in Eq. (|3.17|l below, 
and may be useful as a foundation for the Poisson- Wiseman scheme. Alternatively it may be useful as a check of 
numerical codes that use some regularization scheme. At the core of our analysis is a local expansion of the tail piece 
of the Green's function for linearized perturbations. Such local expansions of Green's functions can be found in the 
literature on quantum field theory in curved spacetime; see Refs. El 

for the scalar case . R ef. for the 

electromagnetic case, and Ref. 57] for the gravitational case. We extend the expansion of Ref. |53 to one higher 
order, and apply the result to compute the quasi-local piece of the gravitational self-force. 

The organization of this paper is as follows. In Section^]we review the formal expression for the gravitational self- 
force obtained by Mino, Sasaki and Tanaka lOj and by Quinn and Wald llj, and define the quasi-local contribution 
to the self-force. In Appendix^ we use symmetry and dimensional arguments to deduce the possible terms that can 
appear in the expansion of the quasi-local contribution, thereby reducing the computation to obtaining one universal 
numerical coefficient at the leading nontrivial order [0(At^)] and four universal numerical coefficients at the next 
higher order [0(At^)]. Section llTTl comDutes those numerical coefficients; the final result is given in Eq. H3.17|l . Some of 
the details of the computation are relegated to Appendices IbI lUl IdI and IeI In Section Hvl we calculate the application 
of the general expression to some interesting cases. Finally, we make some concluding remarks in Section Ivl 

Throughout this paper we use geometrized units in which G = c = 1, and we adopt the sign conventions of Ref. 

lu- 



ll. THE GRAVITATIONAL SELF-FORCE 

A. The Mino-Sasaki-Tanaka-Quinn-Wald formula 

Consider a point particle of mass /i moving on a geodesic x"(t) of a background spacetime (Af, gfj-y), parameterized 
by proper time t. Throughout this paper we assume that the background spacetime satisfies the vacuum Einstein 
equation Rap — 0. The particle will perturb the background geometry to linear order in fi. We denote the linearized 
metric perturbation by haf}, and the more convenient trace-reversed form of this perturbation by 

^Pj = hp^ ~ ^ gp^ hf,u g'"" . (2.1) 
We raise and lower indices with the background metric. We specialize throughout to the Lorentz gauge defined by 

= 0, (2.2) 

where the semicolon denotes a covariant derivative with respect to the background metric gap. In this gauge the 
linearized Einstein field equations take the form of the simple wave equation 

{Dg^ag.p + 2C^aupW = -IfiTrTc,/?, (2.3) 

where □ and C^j,avp are the D'Alembertian and Weyl tensor associated with the background metric gp^, and Tap 
is the linearized stress energy tensor. The wave equation (|2.3(l can be solved using the retarded Green's function 

^ret" ^ 1 which is defined by the equation 

{Ug^ag.p + 2C^a.p)GZ^'^'{x,x') = -g^^a"' gp/ 5\x,x'), (2.4) 
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and by the fact that it has support only when x' is in the causal past of x. Here (/ct" is the parallel displacement 
bivector [sl Is^ l60l| . and b'^(x^x'^ = b^{x — x')/ is a generalized Dirac delta function''. The retarded solution to 
Eq. 1)2. 3|l can be written in terms of the Green's function as 

<tW = Ifivr j d^x' ./'^)GZo.'p'{x, x')T"'P'{x'). (2.5) 
For the point particle source, the stress energy tensor is given by 

/oo 

6^[x,x'{T')]Uo,{T')u0{T')dT', (2.6) 
— oo 

where Ua{T) is the particle's four velocity, and inserting this into Eq. (|2.5() gives 

C;(x) = le^A^ / dr' GZ^.0.[x,x'{t')]u'^' {t')u('' [r'). (2.7) 



-oc 



Now one would expect the particle to move on a geodesic of the total metric gap + i^^^p ~ 9af39''^i''^s /'^- Such 
geodesic motion would be equivalent to motion for which the mass times acceleration with respect to gafs is 

r = M^"''^'^^;%, (2.8) 

where the tensor Pap-yS is given by 

P"^^'' = -iu^it'^u'^M* - g°'^Pu''^u^ + ig^'^u'^w^ + ^u°'g'^^u^ + ^g°'^g'^^. (2.9) 

However, the retarded field V'JfJ and its gradient '^'^^^.^ are divergent on the particle's worldline, so the naive expression 
(|2.8|) for the self force is ill defined. Instead, the correct expression for the self force is given by Eq. H2.8|) with the 
retarded field replaced by the so called tail field ■0*'^'' [13, [lil : 

^^^ail^^^ (2.10) 

Here the angular brackets (. . .) denote the result obtained by averaging over a 2 sphere of some small radius r about 
the particle, in the spatial hypersurface orthogonal to u", and by taking the limit r ^ 0. The tail field V-'Jit''(^) is 
defined by truncating the integral l|2.7|l over the particle's worldline to exclude the contribution from the direct part 
of the Green's function: 

Cii(^) = 16vr/i hm r^^^^ 'dr' GZ.'p'[x, x' {r')] u'^\t')u^' [r'). (2.11) 

«^0+J -oo 

Here Tret (2;) is the value of proper time r at the point where the worldline intersects the past lightcone of the point 
x. The remaining, direct portion of the field is 

^^;ircct(^) = <t(^) -Cii(^) = 16^/^ lim r^'^^^^^ dr' G^:,a'p'[x,x'{r')]u'^' {t')u^' {t'). (2.12) 



Trot (a:)- 



If we now take the gradient of the tail field H2.11|l in order to substitute into the formula H2.1(J|I for the self-force, 
we obtain two terms: a term generated by the action of the gradient operator on the quantity Tret (a;), and a term 
generated by the action of the gradient operator on the retarded Green's function. The first contribution gives an 
expression which has a direction dependent limit on the world line, but which vanishes once the average (. . .) is taken 
[IIII63. The second contribution is continuous on the worldline (so the averaging can be dispensed with), and yields 
for the self-force the expression [Tol [TTI| 

r(T) = 16 VP^'^^Mim r 'dT'Gf^p,^,.s[x,x'{T')]uP\T')u<\T'). (2.13) 



* Note that because we have u sed the sign convention of Misner, Thorne and Wheeler |58ll . our Weyl tensor has the opposite sign to that 
of Mino, Sasaki and Tanaka lid . Note also that our Green's function is defined to be one half that of Ref. llOl . 
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As an aside, we note that Detweiler and Whiting [6l| have introduced an alternative sphtting of the retarded field 
of the form 

i^Z = i^Z, + i^Zn,.r^ (2.14) 

where the singular piece ips^g is a solution of the inhomogeneous wave equation l|2.|-{() which can be computed locally, 
and V'reguiar ^ solution of the corresponding homogeneous wave equation such that the self force is given by the 
expression (|2.8|) with 1/1^^,^ replaced by V'reguiar- "^^^^ alternative formulation also gives rise to the final formula ()2.13fl 
for the self force. 



B. Equation of motion 

The self-force formula H2.13|l discussed above was defined only for a geodesic worldline. Therefore it is necessary 
to supplement the self-force formula with a prescription for computing the motion of a point particle that includes 
the effect of the self force to leading order in the particle's mass /i. The reason that finding such a prescription is 
not entirely trivial is the following llOl the linearized Einstein equation admits solutions only if its source, the 
particle's stress energy tensor (|2.6|l . is conserved. However, the stress energy tensor for a point particle is conserved 
only if the worldline is a geodesic. Therefore it is not straightforward to define the self-force on a non-geodesic 
worldline. 

Quinn and Wald suggested the following method of resolving this difficulty. They define a self-force for a 
non-geodesic worldline by relaxing the Lorcntz gauge condition H2.2|l while retaining the form H2.3|l of the linearized 
Einstein equation. The justification for relaxing the gauge condition is that the associated errors are quadratic in the 
mass /i, while the self-acceleration is linear in /i. Their equation of motion is then 

Ma" = r, (2.15) 

with /" given by Eq. (|2.13(l . modified as described above, and a" is the 4-acceleration with respect to the background 
metric. The equation of motion (|2.15(l is an intcgro-diffcrential equation which is non-local in time. 

Here, we suggest an alternative equation of motion, which gives the same results as Eq. (|2.15f) to leading order in /i. 
First, for any point V in spacetime, and for any unit, future-directed timelike vector u at V, we define the self-force 
vector /" = /"{Vju) to be self force obtained from the prescription H2.13() for the particular geodesic which extends 
into the past from V whose tangent at V is u. The equation of motion for the perturbed worldline (r) is then 

fia^ ^ nx{T),u{T)]. (2.16) 

In other words, the worldline is such that its acceleration at any point is the acceleration that is obtained from 
the integral (|2.13|) for the geodesic which is tangent to the worldline at that point. The justification for using the 
instantaneously-tangential geodesic rather than the true worldline is essentially the reduction of order argument 
discussed in Refs. [ill|62|. The equation of motion (|2.1(j|) is a second order differential equation which is local in time. 



C. Definition of tlie quasi-local piece of the self- force 

In this paper we evaluate not the entire integral over the worldline in the self- force formula H2.13|l . but instead that 
portion of that integral near the particle with proper times r' in the range r — At < t' < t. Specifically, we define 

/^L(r,Ar) = 16VF"'5'^Mim r ' dr' G^p^^f,,^,.^s[^,x' {T')]uf'' {r'^' (r'). (2.17) 

O+J T-Ar 

Here the subscript "QL" denotes "quasi-local" . By comparing with Eq. H2.13|l we see that the entire self-force is 
obtained from /ql(t. At) in the limit where At 00. We will derive an approximate power series expansion of 
/ql(t. At) that is valid in the limit At — s- 0. 

Within a sufficiently small neighborhood of the point a;" = x°'{t), the retarded Green's function can be written in 
the Hadamard form 

C^r'^'ix^x') ^ ^&^x{[/^''"''5'(a:,^') S[a{x,x')] - V^'"^'^' (x^x') e[-aix,x')]}. (2.18) 
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Here T,{x, x') is an arbitrary function which is positive for x in the causal future of x' and negative otherwise, 8[-] is 
the Heaviside step function, a{x,x') is Synge's world function 0,|63|, is the ordinary Dirac delta distribution, 
and U^"'°'''^' {x,x') and V^"'°'' 1^' {x,x' ) are smooth functions. The part of the Green's function proportional to 5{a) is 
called the direct part, and the part of the Green's function proportional to the function V^'^'^ ^ as the tail part. For 
sufficiently small Ar, the portion of the worldline between x" (r — Ar) and (r) that arises in Eq. (|2.17|l will lie 
inside the neighborhood where the Hadamard expression H2.18|l is valid. Inserting this expression into the formula 
(|TT7|) for gives 



/^L(T,Ar) = -4M2p 



Vp^p,y.s[x,x'{T')y {t')u^' {T')dT' . 



(2.19) 



t-At 



Note that the direct part of the Green's function does not contribute to the expression (|2.19(l , because of the limiting 
process involving e in Eq. (|2.17(l . That limiting process is unnecessary for the tail contribution which gives an integrand 
that is finite at r' = r; hence there is no limiting process in the final result H2.19|l . 

The integral H2.19() can be approximately evaluated for small Ar using a covariant local expansion of Va-y/3'-y'(a;, x') 
in a neighborhood of x' — x. The result is of the form 



/Qi„(r, Ar) = /(")(r) + fi'\T)AT + fi'\T)AT' + .eHT)AT^ + 0{At'). 



(2.20) 



Because this expansion arises from a covariant Taylor series, the coefficients fa^ (r) are purely local geometric quan- 
tities evaluated at the present position x"(r) of the particle. In Appendix 1X1 we derive the geometric content of these 
coefficients using simple counting and dimensional arguments. We find that 



f(°) = 

J a 

f(i) = 



(2.21) 
(2.22) 
(2.23) 
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+C3 



V A 



V ;A 

<5/3 



(2.24) 



where co,ci, . . . ,C4 are dimensionless numerical coefficients that are as yet undetermined. The calculation of these 
coefficients proves to be the most difficult part of determining the expansion of /q^:,, and is the topic of the next 
section. 



III. EXPANSION OF THE SELF-FORCE 



In this section we use a local covariant expansion of the retarded Green's function to compute the numerical 
coefficients appearing in the expansion H2.20|l of the self force. We use the formalism of bitensors developed by 
DeWitt and Brehme see Poisson |l3| for a recent detailed review of this formalism. A fundamental role in this 
formalism is played by Synge's world function g(x,x'^ and its derivative (7;Q,(a;, a;') (see Appendix Ib|| . The local 
covariant expansion of any bitensor T{x, x') takes the form 

T(x, x') = V -tl^-<^-{x)G,^, (x, x')... a,^^{x, x'), (3.1) 
■^-^ n! 

n=0 

where the coefficients t"^ - "" are local tensors at x. For bookkeeping purposes we define = \a\, then it follows from 
Eq. ljB5|l below that the nth term in Eq. (|3.1|l scales as s". We shall use s as an expansion parameter throughout our 
computations. 

As a foundation for the expansion of the retarded Green's function, we compute in Appendix ^ local covariant 
expansions of a number of fundamental bitensors, including the second derivative a-^ai3 of the world function, various 
covariant derivatives of the parallel displacement bivector g^" , and the Van-Vleck Morette determinant A, to order 
O(s^) beyond the leading order. These expansions were originally computed to O(s^) by Christensen [sol l60l|. and 
extended to O(s^) by Brown and Ottewill 54, "s^. Our results agree with those of Brown and Ottewill, except for 
one case where we correct their result [Eq. (|B33|) below]. 
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In Appendix [HI we compute the expansion to order Ofsf) of tail portion Vaf3a'i3' of the retarded Green's function, 
extending previous work of Allen, Folacci and Ottewill |57j who computed the expansion to order 0{s^). The result 



is of the form 



(3.2) 



cf. Eqs. (|C9|I and l|C14(l below. Here the various expansion coefficients w° ^(a;) and w„ ^(a;) are given in Eqs. I|C15|I 
- (tni^ and (in22ll - (in23ll below. 

We now turn to evaluation of the integrand in the expression H2.19|l for the quasi-local piece of the self force. First, 
we note that the four velocity is parallel transported along the worldline, so we can make the replacement 



We can rewrite the first factor on the right hand side as 



(3.3) 



(3.4) 



where we have defined Vap-yS = g^" g^^ Vapa' p' and we have used the definition IjBlSp of the tensor Qap-y- Using the 



expansions (|3.2|l . (|B15|) . and (|B25|) we now obtain 



where 



^0 A ^0 ^1 -^0 



(3.5) 

(3.6) 
(3.7) 

(3.8) 



It is understood that the right hand sides of Eqs. 1)3. 6|) - H3.8|l are to be symmetrized on the index pair {^xv), and 
on the index pair (pr) if present. We also note that a'°'{x,x') is proportional to the tangent to the geodesic joining 
X = x°'{t) to x' = x" (r'), i.e., the four velocity u". It follows from the normalization condition (|B5|I that 



a' {x,x ) = [T ~ T )u 



(3.9) 



where is the four velocity at x"(t). Finally, substituting the formulae H3.3|) . (|3.5|l and l|3.9ll into the expression 
(|TT^ gives 



/^Jr, At) = -4^2 P^^^' u^u^ / dr' {v^,^., + V^^^Jup^r r') + V^.^^,/^ UpuAr - r'f +0[{t- r'f] } . 

(3.10) 



Evaluating the integral over t' gives an expansion of /q^ of the form H2.20|l . where the coefficients are 

y(0)a ^ 



(3.11) 
(3.12) 

(3.13) 
(3.14) 



The various coefficients Vq...^ are obtained from the formulae (|3.6() - (|3.8() and are tabulated in Appendix IdI We now 
substitute those expressions into Eqs. H3.11|l - H3.14|l and use the definition (|2.9|l of the projection tensor p"/57<5_ After 
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a considerable amount of algebra we obtain coefRcients fa^ of the form H2.21|l - (|2.23|l . as expected. The numerical 
values of the coefRcients are 



Co 



-1 



and 



Cl = -, C2 
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20' '^^3' 



19 
'60' 



(3.15) 



(3.16) 



Our final expression for the quasi-local piece of the self-force is therefore 



19 
60 



>Ar3 + 0(Ar4 



(3.17) 



IV. SOME SPECIAL CASES 



Our expression H3.17|l for the quasi-local contribution to the self-force is quite general, applying to a massive particle 
with any 4- velocity in any vacuum background spacetime. There are, however, cases which are of more intrinsic interest 
than others. In particular, cases in which the background is a black hole are of interest for developing templates for 
LISA. We now examine the form that /q^ takes for several particle four-velocities in a Schwarzschild background and 
for two simple four-velocities in a Kerr background. 

We begin with the Schwarzschild background. In standard Schwarzschild coordinates, the line element is 



+ r^de^ 



9 • 9 n 1 i9 



(4.1) 



where m is the mass of the Schwarzschild black hole. Consider a general particle four velocity. Because of the spherical 
symmetry of the background, we can always arrange that at any instant, the particle as well as the spatial projection 
of the particle's 4- velocity be in the equatorial plane. Doing so, we have that 9 = Tr/2 and = for a general particle 
4-velocity. This leave three non-vanishing components for u", however, the normalization condition u°'ua = —1 can 
be used to eliminate one (we have chosen to eliminate u*). Thus, in this background, the velocity of of the particle is 
completely specified in general by a choice of = dr/dr = f and u"^ = dip/dr = 0. 
Putting this general 4-velocity and the Schwarzschild metric into (|3.17(l yields 



m 



{r — 2m) 



(^r + (r — 2m) + 1^ ^ 



9 02 ^2 02 ^2 _,_ 1^ + ^ f (^64 + 150 r'^ + l) At^ + ©(Ar"* 
- 9r* {2 it? r2 -I- 1^ 



■r At^ 



20 ' -I- 80 TO (/j"* r"* + 4 m - 40 0"* - 2 r f2 



(4.2) 
(4.3) 
(4.4) 



- 300 f2 0^ - 2 r - 128 0^ - 31 0^ \ Ar^ -\- 0{At'') 



„3 -2 12 



12 „3 



/ql 

f<t> 



9 r (j)^ + 1) (2 02 r2 + 1^ Ar2 -|- ^ (388 ij? + 300 0^ + 99^ r At^ + 0(At*) 



(4.5) 
(4.6) 



One expects, of course, the vanishing of the 9 component by symmetry arguments. However, we notice several other 
features. For purely radial motion (0 = 0), the quasi-local part of the self-force vanishes to O(Ar^), and the /q^ 
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component vanishes to 0{At'^) - indeed, by symmetry arguments it must vanish to all orders. Thus, the quasi-local 
part of the self-force is directed toward the black hole to order 0{At'^) in this case. For purely tangential motion 
(r = 0), /ql vanishes to O(At^). Further, there is no Ar^ contribution to /ql or /q^ in this case. On the other 

hand, /g^ does not vanish even if for a static particle, where ^ = r = 0. 

A particularly interesting case is that of a particle following a circular geodesic around the black hole. In this case, 



which gives a quasi-local self-force contribution of 



.3 



ft _n„2^ / ' jr-m) , ^^^_4^ 

r6 ^ r-2,m (r-3mf 

f ^ 1 (35 - 19 mr - 2 r^) Ar^ + O(Ar^) (4.9) 

/^L = (4.10) 



/I = - f (. - 2H A.^ + 0(A.^) (4.11) 



Interestingly, we can assign physical meanings to two of these components. Assuming that the self- force causes an 

adiabatic deviation from the background geodesic, we can define an energy for the particle of E = ut and an angular 
momentum of L = u^. Thus, dE/dr = du^ /dr = f*^ / jJL and dL/dr = dy'^/dr = f^/n- In other words, /ql/a* a^^d /ql 
can respectively be interpreted as the energy and angular momentum radiated by the particle due to the quasi-local 
part of the self-force. 

Of course, we expect astrophysical black holes to provide a Kerr background spacctime, so this background is more 
astrophysically relevant than Schwarzschild. Unfortunately, as is often the case, calculations in the Kerr background 
lead to longer and less managable expressions. Using symbolic algebra programs, it is straightforward to calculate, 
for example, for the case of general equatorial motion in Kerr, and we have done so using MAPLE. The expressions, 
however, are so unwieldy that we choose not to display them here. Rather, we illustrate with the two simplest motions 
in a Kerr background, stationary with respect to an observer at rest at infinity and co-rotating with the black hole. 
We restrict the particle to be in the equatorial plane in both cases for convenience. 

We start with the metric in Boyer-Lindquist coordinates. 



, 2 / ^-a^sin^6'\ , , ■ 2 r. f r'^ + ^'^ - ^\ , , 

ds^ = - dt^ - 2osin2 9 dtdcp - 



where 



(r^ -l-a^)^ - a^Asin^( 



sin^ 9d(h^ + ^dr^ + ^d9'\ 



(4.12) 



S =r^ -l-a^cos^fil, (4.13) 

A=r'^ -2mr. (4.14) 

As usual, a is the spin parameter for the Kerr black hole and rn is its mass. Wc will only bo considering the 
special case where the particle is located at = 7r/2, and only particles with = =0. Again, we have the 
normalization condition u^Ua = — 1 which we use to further eliminate u*. We therefore have only one specifiable 

4-vclocity component, (f>. 

We begin with the case where ^ = 0. This is the case of a particle which is momentarily stationary with respect to 
an observer at spatial infinity. In that case, the quasi-local part of the self-force becomes 

(2\ 3/2 
r (r - 2 m) (r^ - 2mr + a^) 

^QL=-M -l^--^^ ^-—-4 ^Ar^ + 0(AA (4.15) 



r 



(r^ — 2 mr + a^) (r — 2 m) 



/^i^ = - A ^2 ^ 2 mr + a^) ^ ^ - 4 r^m + + 4 mV^) At^ + 0{At''), (4.16) 

f=0, (4.17) 

3/2 



(j.2 _ 2 mr + 2 a2) (r {r - 2 m) (r^ - 2 mr + a^) 
{r'^ -2mr + a^f {r-2mf 



ft. =-^^' ;::::;3 — — ^^-^ — ^ + o(Ar^). (4.i8) 



10 



It is interesting to note that there is a self-force at order Ar^ in this case, unHke the case of a static particle in 
Schwarzschild. 

One might argue that this is not unexpected, since in Kerr the particle is "moving" with respect to the rotating 
background. A fairer comparison, therefore, might be with a particle that is co-rotating in the Kerr background, i.e. 
one for which = 0. However, in this case wc find the quasi-local part of the self-force is 

--m'^-^^^^^^^^^^-^^^ f22aWr2-27aVm-71aVm-44a6^m (4.19) 
^ 10 (r3 + 2ma2 + a^r) V 

+ r** + 10 r^a^ + 27 aV^ + 28 a^r^ + 10 Ar^ + 0(At'') 

f^^^ - 9 ^ / (4.20) 

r (r^ — 2 mr + a^) {r^ + a^r + 2 ma^) 

X {r{r^ + a\ + 2ma^) {r^ ~ 2mr + a^)f^^ Ar^ + ©(Ar^) 

Thus, in both analogues of the static particle in a Schwarzschild background, the rotation of the Kerr background 
induces a radiation reaction force at order At^. 

We end this section with a warning. While the expressions derived here might be useful for comparisons, they do 
not, in general, have any intrinsic meaning (the case of the circular geodesic in Schwarzschild being somewhat of an 
exception) . One obvious reason for this is that we have found only a part of the self- force which reflects a very limited 
part of the particle's world-line. However, there is a more subtle limitation as well. The self-force is a gauge dependent 
quantity, depending on the perturbation gauge chosen. Our expressions, which are composed of quantities that are 
gauge invariants of the background, are nonetheless tied to the Lorentz perturbation gauge. Thus, these expressions 
may, in general, only be legitimately compared to other expressions for the self-force in the Lorentz gauge. 



V. CONCLUSION 



In this paper, we have discussed a novel approach to calculating the self-force experienced by a massive particle 
moving on a geodesic in a curved background. In this approach, proposed by Poisson and Wiseman, one does not 
regularize the retarded Green's functions nor the retarded field. Rather, one explicitly calculates the tail part of the 
retarded Green's function in the normal neighborhood of its current position. From this, one can calculate the "quasi- 
local" part of the self-force, which arises from some finite portion of the particle's worldline within this neighborhood, 
of duration Ar. The rest of the self-force can then be obtained using the full retarded Green's function (without 
need of regularization) , since it is well behaved (and gives the correct contribution) for the portion of the particle's 
world-line beyond Ar to the past. 

We have also carried out the first step in this procedure, the calculation of the first two non-vanishing terms of 
an expansion of the quasi-local part of the self-force in Ar. Our expression has some remarkable properties. First, 
it is quite general, in that it does not rely on prior specification of the particle motion nor on prior specification 
of the background geometry (we have restricted our attention to vacuum backgrounds in this paper, but this was 
for convenience, and exactly the same procedure can be used for calculating the quasi-local part of the self-force in 
backgrounds with matter). Second, we express the quasi-local part of the self- force in terms of quantities that are 
purely local to the particle. Thus, our expression does not require a detailed understanding of the past history of the 
particle - it is only a function of the current position and velocity of the particle. This can hardly be surprising - 
within a normal neighborhood, there is a unique geodesic specified by any 4- vector at a given point, which the particle 
is assumed to be travelling along. 

Much work remains in order to determine even if the Poisson- Wiseman prescription is feasible, let alone to calculate 
a complete self-force using it. While the approach here is quite general, it is likely that the retarded Green's function 
from which the rest of the self-force is calculated will be most easily obtained in most scenarios of interest by a 
mode-sum expansion. It remains to be seen whether it is technically feasible to calculate the quasi-local part of the 
self-force to sufficient order that it has sufficient precision at a distance from the particle at which the mode-sum 
converges well. Work is currently under way to begin addressing such questions '6(f|. 

Some hope arises from the work of Anderson and Hu(67|, who have shown how to calculate the tail part of the 
retarded Green's function for a scalar particle in Schwarzschild using the Hadamard-WKB approximation. This ap- 
proach, which should work for spin-1 (electromagnetic) and spin-2 (gravitational) fields on a Schwarzschild background, 
allows one to calculate to a much higher order in Ar for a fixed amount of effort. Furthermore, it might be extendible 
to Kerr backgrounds. In any case, we echo their sentiment, that this approach warrants further investigation. 
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APPENDIX A: GEOMETRIC CONTENT OF THE EXPANSION COEFFICIENTS 

In this appendix, we use symmetry and dimensional arguments to deduce the forms of the coefficients f!^^\ fa\ 
fa'^ and fa"^ appearing in the expansion (|2.20|l of the quasi-local piece of the self force, up to some unknown numerical 
coefficients. We start by noting that the self force is proportional to the square fj? of the particle's mass, from Eq. 
(|2.19|) . Factoring out this factor of /i^ , we can write the expansion (|2.2U|) as 

/Qi„(r, At) = ^? [/(°) + /(I) Ar + fl^'^Ar' + fi^^Ar' + 0{Ar')] . (Al) 

Here the coefficients fa^ — fa^/^P satisfy the following key properties: 

• They are independent of the particle mass fi. 

• They must be constructed as polynomial expressions in the following tensors at x"(t): the metric gap, the 
four velocity m", and the Weyl tensor and its various symmetrized derivatives Cap^s, Caf}-yS;e, Cai3jS-{ep)j 6tc. 
The Ricci tensor does not appear since we are assuming a vacuum background, and the four acceleration of 
the curve a" does not appear since we are assuming a geodesic curve at zeroth order. The Levi-Civita tensor 
eQ/37(5 cannot appear, since the expression H2.19|l for the self force is invariant under the parity transformation 

• Since the force is dimensionless (in geometric units in which G — c ^ 1), and At has dimensions of length, each 
coefficient fa^ has dimension (length) 

• Each coefficient fa'' must be orthogonal to the four velocity m", since the total force (|2.19|l has this property. 

We now apply these properties to deduce the most general allowed forms of the various coefficients, which are given 
in Eqs. ifT^ - above. 



1. The zeroth order coefficient /q"' 

The coefficient fa'^ has dimension (length)"^, and hence must be linear in Cap-yS which also has dimension 
(length)"^. None of the derivatives of the Weyl tensor can appear. However, there is no non- vanishing vector 
that can be constructed out of Cap-yS, 9ap and u" that is orthogonal to u". One needs to contract at least three 
of the indices on the Weyl tensor with the metric or the four velocity. However the Weyl tensor is traceless on all 
pairs of indices, so one cannot use the metric to contract any pair of indices. Also the antisymmetry properties 
(^{apj-yS = C!ap{-yS) = of the Wcyl tensor mean that one cannot contract with three factors of four velocity. Hence 

the coefficient fa'' must vanish, cf. Eq. (|2.21l) above. A version of this argument was first given by Ori, and was used 
to deduce the fact that the formula for the gravitational self force could contain only the tail term and could not 
contain any local terms, unlike the scalar and electromagnetic self- force expressions 0|. 

An alternative, simpler version of the argument can be obtained by considering the independent, electric and 
magnetic components of the Weyl tensor. Introduce an orthonormal basis e? = u" and e?, 1 < j < 3. Then we define 

£i]^CQm^e-e^^elCap,s (A2) 
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and 

% = -\^iyGu(^^ - \e.p,,C^\,u'^4u^e]. (A3) 

In vacuum these are symmetric, traceless tensors, and the Weyl tensor can be expressed in terms of these tensors via 
the formulae 

= -m^ii (A4) 

and 

^Ijkl — ~^ijp^kiq^Pg (A^) 

together with Eq. ljA2p . The coefficient /i"-* must be a three- vector that is linear in fj--, or linear in B-^j with one 
factor of e-^jf, by parity arguments, and there is no such vector since yB^j-j — 0. 



2. The first order coefficient /i^' 

The coefficient fa"^ has dimension (length)"'^ and so must be linear in Cap-fS-.e- Again, it is impossible to form a 
non- vanishing 4- vector that is orthogonal to by contracting Cap-ys-e with the metric and/or four velocity. The 
additional derivative index e cannot be contracted with any of the indices on the Weyl tensor, since by the Bianchi 
identity the Weyl tensor is divergence free on all its indices in vacuum: 

The e index can be contracted with the four velocity to form Cap-yS-eu'^ , but then one is faced with the same problem 
as above of contracting three of the four remaining free indices to obtain a vector. It follows that fa ^ must vanish, 
c/. Eq. above. 

One can also phrase this argument in terms of the electric and magnetic components of the Weyl tensor as before. 
The components of CagyS-e can be represented as the "time derivatives" and "spatial derivatives" of £-■-- and B~-~-: 

£^ EE C^pys,,u'^e?u^e'.u' (A7) 
^ Cc.p,5;eu''elu^e\el (A8) 

% ^ -\e^p,8C-",^.,^u-eP^u^e^^ur (A9) 
^5],^ ^ -\eo.0,8C^',,,„u'^e^u^e^^el. (AlO) 
These quantities are not all independent but obey constraints that follow from the Bianchi identity: 

- 0, (All) 
B^- ^ 0, (A12) 

% = -^mHj^ (A13) 
% = ^lid^m- (A14) 

Since the time derivatives can be obtained from the spatial derivatives, we can without loss of generality restrict 
attention to the spatial derivatives. Thus, the coefficient /i^^ must depend linearly on £^ f, and/or Sj- f,, and can 
depend in addition only on and on e^-.^. Consider first the derivative This quantity cannot be contracted 

with e^-j, by parity arguments, and it is easy to see that one cannot obtain a nonvanishing spatial vector. The only 

candidate vectors are the divergence E-^j j,5^^ which vanishes by Eq. ljAll|l . and the contraction £-q ^,(5*^ , which vanishes 
by the traceless property of E^. Next, consider the derivative f,. By parity, this quantity must be accompanied 
by one factor of e^-.^ (or an odd number of factors of it). Any such tensor will have an even number of indices, and 
consequently it is impossible to obtain by contraction a vector. 
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3. The second order coefficient fa 

The coefficient fa'' has dimension (length)^*, and so must be either quadratic in Cap-yS or finear in Cap^s-Xea)- 
Consider first the second derivative term CapjS-^ea) ■ Let us first analyze the term Cap-yS-ea- without the symmetrization. 
The six indices on this tensor must be reduced to one index by contractions with the metric and/or the four velocity. 
As before, the problem is getting rid of at least three of the indices a, f3, 7 and S. The only new type of contraction 
that is available is to contract one of these indices with the second derivative index a [contractions with the first 
derivative index e vanish by Eq. (|A6p ]. However, such a contraction can be re-expressed as a product of two Weyl 
tensors by commuting the indices. The same argument also applies to the symmetrized derivative Caf3yS;{sc7): since it 
is a linear combination of un-symmetrized derivatives. Thus, it is sufficient to consider expressions that are quadratic 
in the Weyl tensor ^. 

Consider therefore expressions that are quadratic in the Weyl tensor, or linear in Cap-ysCeapX- We can classify such 
expressions in terms of the number of contractions between indices on the first Weyl tensor and indices on the second 
Weyl tensor. The cases of zero, four and one contractions are easy to dispense with. For example, in the case of one 
contraction, there is no vector orthogonal to that can be obtained by contracting Cap-ysC^^p with g^p and u" 
(without further contractions between the two Weyl tensors). In the case of three contractions, there are two different 
tensors that one can construct, namely 

Vap = CaySeCp^'' (A15) 

and 

Wap = Ca^seCp''''. (A16) 

However it follows from Ca[py5\ — that Vap — 2Wap, so there is only one independent tensor obtainable from three 
contractions. Furthermore using the formulae (|A2|I . ljA4|l and (|A5|I one can show that Wap is proportional to the 
metric: 



Wap^-Wy'^gap- (A17) 



Hence there is no non-vanishing vector that can be formed using three contractions. 

Consider next the case of two contractions. There are three different four index tensors that can be obtained with 
two contractions, namely 

Xapea — C'ap^sC^^"^^ , (Al8) 
Ya^^p — Cap^gC^^p^, (A19) 

and 

Zayecr = CaPysC^^'^^ . (A20) 

The tensor Xapea does not yield any candidate expressions, since one needs to contract it with three factors of the 
four velocity, and it is antisymmetric on the index pairs (a/3) and (sa). Similarly the tensor Za-yea is antisymmetric 
on the index pairs (07) and (ecr), and so does not yield any candidate expressions. The tensor Ya-yep can be used to 
construct the non-vanishing 4-vector Ya-yepu'u^'^'' ■ Projecting this vector orthogonal to the four velocity yields the 
possible term 

{5i + Uau'')Cpys,C,\'u'u''uP. (A21) 



In fact, the only nontrivial candidate expression one can construct from the symmetrized second derivative vanishes. This expression is 

CaMS;(ea)9'''uPu-'u'^ = -\ [c .^cT C ^0-," + C ^[ C ^.y" + C pC ^^f] g^'uf^ulu^, 

where we have used Eq. iA6i . Here the last term in the square brackets vanishes since the the first factor is symmetric in the index pair 
(/!>(t) by virtue of being contracted with u^u"' , and the second factor is antisymmetric in (per). The remaining two terms in the square 
brackets cancel against each other. 



14 



This quantity is nonvanishing in general and and so can appear in the expression for fa , cf. Eq. H2.23|l above. It is 
the only term that arises at this order in At. Using the formulae ljA2|l and ljA4|l . we can express this term in terms 
of the electric and magnetic components of the Weyl tensor as 



(A22) 



"(2) 

This completes the derivation of the most general allowed form of /A . 

The analysis of expressions that are quadratic in the Weyl tensor can be rephrased more simply in terms of the 
tensors fj; and B^. One needs a spatial vector that is bilinear in and/or B^epqf, since by parity each factor of 
B-fj must be accompanied by a factor of epqr- It is easy to see that the only non- vanishing candidate expression is the 
product (|A22|I . 



4. The third order coefficient fa 



(3) 



The coefficient fa^ has dimension (length)"^, and so must be either bilinear in CapjS and Ca/i^yS-e, or else or linear 
in the symmetrized third derivative Caf3jS;{e<7p)- The argument used in the first paragraph of Sec. lA 3l above also 
applies here and shows that any expression constructed from the third derivative can be expressed as a product of 
the Weyl tensor and a first derivative of the Weyl tensor. Therefore it is sufficient to consider such products. 

For analyzing these products, a fully covariant analysis would be very complex, so we use the simpler formalism of 
the electric and magnetic components. We need to construct a vector that is a contraction of a product of one of the 
tensors 



together with one of the tensors 



B-fjipgr 



(A23) 



(A24) 



As before we can neglect the time derivatives because of Eqs. (|A13I) and IIA14|) . It follows from the identities HA11|) 
- ljA14|l and the fact that all of the tensors (|A23p and (jA24p are symmetric and tracefree on the index pair {ij) that 
there are only four possible nonvanishing vectors that one can construct, namely^ 



'^jk'^ij.k' 



B-~B-~ - 

'-^jk'-^jkA^ 



and 



B-~B-~ ~ 



The corresponding covariant expressions can be obtained from Eqs. HA2() and HA3|I and are, respectively, 



(5f + Uauf')C\^s.C//.^u^u'u'u' 



5f + uau^)u''u^ 
and 



S ;/3 



(A25) 

(A26) 
(A27) 

(A28) 
(A29) 



cf Eq. above. 



APPENDIX B: LOCAL EXPANSIONS OF SOME FUNDAMENTAL BITENSORS 

In this appendix we compute local covariant expansions of a number of fundamental bitensors, including the second 
derivative a-a/s of the world function, various covariant derivatives of the parallel displacement bivector g^" , and 



One can also construct from the time derivatives the possible terms ^iji^t ^jl^ii^ and Hjk^ji^ik' However these terms can be expressed 
in terms of the four quantities iA25i using the formulae IA13i and 



ijk J I Ik 



jk jk,t 



jk ij,k' 



e^-.:B-f£:f 
ijk J I Ik 
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the Van-Vleck Morette determinant A, to order O(s^) beyond the leading order. These expansions were originally 
computed to 0{s'^) by Christensen US 113, and extended to O(s^) by Brown and Ottewill 54, 55^. Our results agree 
with those of Brown and Ottewill, except for one case where we correct their result [Eq. (|B33|I below]. All of the 
formulae in this appendix are valid in any number of spacetime dimensions, except in specific cases, noted below, 
which make use of the identity (|B4|I which is specific to four dimensions. For a detailed and pedagogic review of the 
computational methods used here see Ref. fi^ - 

We first note a number of useful identities that follow from the Bianchi identity and from properties of the Riemann 
tensor: 

RapisRe^^'^ -.^j. = ^^aPjsRe'^^^ -p, (B2) 
Ra'^'^^R/lfj.jv.X = —-^Ra^'^^RlStivXil (B3) 

Cafii'xCij^'^^ — -gapC^iyXpC^'"^'' . (B4) 

Next, we derive the coincidence limits a;' ^ a; of symmetrized derivatives of the world function a. Following Ref. |^, 
these can be obtained by repeatedly differentiating the identity 

<^ = Ig'^^^-.o.^-.P (B5) 

and by taking the coincidence limit. The results are 

lim a.ai3 = gap, (B6) 



lim a-a/B-y = 0, (B7) 
lim a-aPfS = --r [RajfIS + RaSf3j] , (B8) 

X —^x o 

lim a-al3iSe = ^7 [RajISSie + RaSlSf.e + RaSf3e;j + RaeflSif + Raepi\5 + Ra-yfie-s] , (B9) 
X — ^x 4 

12 8 

lim (T;Q/3(75e^) = — -r-i?a7/35;eM — T^Ra-rSxRlSen ' (BIO) 

and 

^}^^(^;a(3{'ySefiiy) = -^Ra-tl3S,epiu -^^a-ySxRpeJ^ -V (BH) 

In Eqs. IjBlOp - ljBll|l . the right hand sides are to be symmetrized over the index pair {a(3) and over as many of the 
indices {■ySenu) as are present. 

Next, consider any smooth bitensor T{x, x') (we suppress tensor indices on T). We can expand this bitensor about 
X as a covariant Taylor expansion 

OO - 

T{x, x') = V {x)o,o., {x, x')... a;«„ {x, x'), (B12) 

n=0 

where the coefficients - are local tensors at x. By repeatedly differentiating Eq. (|B12p and taking the coincidence 
limit we find for n = that 

to{x) = lim T{x,x'), (B13) 
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together with the recursion relation for n > 1 

inai...Q„ (2;) — 2^;(qi...q„) (2^1 2; ) ^ ^ ( r ) i°'^■■■'^r^,ar+l...an)i■^)■ 



n-l 



r=0 



(B14) 



We now apply the recursion relation (|ljl4f) to the bitensor a-apix, x'), using the formulae (|lj6f) - HB11|I . The result is 
the expansion 



1 

120" 



1 



It follows from Eq. I)B15|I that in vacuum 

Da = n + 0{s^), 
where n is the number of spacetime dimensions, and that 

□ (cr^") === 0(s3). 

A similar computation can be carried out for the tensor Qap^y defined by 
By starting with the identity ^ 

repeatedly differentiating and taking the coincidence limit one obtains 

lim ga'a-p = 0, 



^;7^;<5^;e^;M^;'' + 0(56)^ (B15) 

(B16) 
(B17) 
(B18) 
(B19) 
(B20) 



lim ga'a-ii-t g"" p — —-^RapPi, 



(B21) 



1™ 9a'a;3('y5) 9^ p — 



(B22) 



lim ga'a;!3(^Se) g p — ^Rapf3^\5e ^R-ap^^pR^ 



and 



y a' __R _R R _ _ ?? R 



(B23) 



(B24) 



In Eqs. (|B20() - HB24|I . the right hand sides are understood to be symmetrized over as many of the indices {'ySen) as 
are present. Combining the coincidence limits (jB20|) - (jB24|) with the formulae ljB12|l - ljB14|) for Taylor expansions 
together with the definition IjBlSp one obtains 



QaPi — -^Ra0fS<^'^ ~ -^Rap-tS:,e'^'^ (^''^ + ^ [RaP-fS-.ep + Raf}5pR~)ep''\ Cr'^a''^a'^ 



1 

120 



7 p ^ p 

Ral3^S;epu + '^^RapSp;uR-yep ^" '^RaPSpR^ep, -.v 



a'-^CT'^a'^a''' + O(s^). (B25) 



Next, we derive an expansion for the quantity g^" ga'B-^'^ by using the identity 



(B26) 
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together with the expansion ljB25|) . The resuh is 

"1 7 2 
_ri P V I _r< ry P V , _n ri P V 

where we have speciahzed to the vacuum case Raf} — 0. 
Next, following Ref. Q we define the tensor 

Dap = -g" ^\a'l3 
which is related to the Van Vleck-Morette determinant [g^ A by 

A ^ detL>/. 

Using the identity 

Dap ~ QajP Cr''^ + Cr-ap 

together with the expansions (|B15|) and (|B25I) we obtain the expansion 

1 



a'"'a-^a-' + Ois'^), (B27) 

(B28) 
(B29) 
(B30) 



Dap — 9aP + QRajPSCr''^a'^ — —RajPS-e'^'"' O^'^ a''^ 



40 



1 



1 



1 



^a-yop-^^ 

The determinant of ^ can be calculated via (|E4|) . Taking the square root of the determinant gives 



^^QRa-fPS-Efii' + ggRa-jSp^uRp^ff + ^2Q^°'^^P^0'^J' -y 



(B31) 



7 



12 
1 



1 



1 



80 



RaPijS + -Z^RpaP^R -yS + l^RapRjS 



1 



360 ^"'^''^''^ ^ 288 "^"'^^'^^''''^ ~'~ ^QQ^P°'^'^^''i^'^'<^ 



ggQ--p^^v^-- 288 



■,a -.(S ■'J -S 



(B32) 



Acting on this expression with the wave operator and using the identities ljBl|) - (|B3|I together with Eq. ljB15|l we 

obtain 



□Ai/2 ^ \r 
6 



— - -^R;ap + -^RRaP " ^^a7^''/3 + ^^07/3*-^'''^ + '^Ro^lSeRp''^^ 



^\a-p 



— R R — R R P ^ _ R 'P _ D R P 

^^rLapiti^ 180^" " '"^ 120 ■'"'^ 180 ^'''" '^ 



•-P 7 



g^RapuxRp^ ;7 i44-^-^"'3;'>' + 360^'"'''^ 
Specializing to the vacuum case Rap — yields 



In the special case of 4 spacetime dimensions this can be further simplified using the identity ljB4p to give 
□Ai/2 = ^C.^.aC^^'^' 5a/3 a^"(T^'5 - -^C^^xp^.C^-^Pgapa^'-a^^cj-^ + ©(s^)- 



(B33) 



(B34) 



(B35) 



APPENDIX C: COVARIANT EXPANSION OF THE RETARDED GREEN'S FUNCTION 

In this appendix we derive the local covariant expansion of the retarded Green's function Gj^'t" ^ (x, x') which is 
defined by the differential equation 



{^gt^agvP + 2C^a,.p)G'^^^'^ ^ (X, X) 



gia"^' gp/ + ugapg"'^' 



S*{x,x'). 



(CI) 
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Here we have introduced a real parameter k to facilitate comparison with the work of Allen, Folacci and Ottewill 
(AFO) [s^l, who analyzed the case k = —1/2. For this paper we are interested in the case k — 0, cf. Eq. (|2.4|l above. 
Note that these two Green's functions, the cases k = and k = —1/2, are related to each other by a trace reversal on 
the index pair (a'/3'). Throughout this appendix we specialize to 4 spacetime dimensions. 

We use the standard method explained by Hadamard 5], DeWitt and Brchme^Sj, and AFO js^l in the scalar, vector 
and tensor cases respectively, and we extend the expansions of AFO to one higher order. We assume for the Feynman 
Green's function the expression 

G^" ^3 (x, x') = A ( + y'^"'" ^ In \a + ie] + VF'"'" ^ (x, x') \ (C2) 

for some bitensors f/^'''"'/'', yM^"'/3' and W>'''°''^\ where we have introduced a regularization parameter e to give 
the appropriate singularity structure. The expression H2.18|l for the retarded Green's function can be obtained by 
taking the negative of the imaginary part of the Feynman Green's function l|C2|) . and by multiplying by the function 
e[E(x,a;')] defined in Sec.lnTl!^. 

Substituting the real part of the Green's function ljC2p into the homogeneous version of the differential equation 
(jGl|) . and equating to zero the coefficients of I/ct'^, Incr, and the remainder gives the three equations 

f;MJ.a'/3';7^.^ _ 1 JJP-^a'/J' (J^ A)''^a;^ = 0, (C3) 



^a/S^'^V^.a'^' = 0, (C4) 

and 

2)0/3^" C/^.a'/3' + 2V;,/3a'/3' + 2K.^„,^/V;^ - V^(i^, p, ^)'^ U ,^ + d 2)„/3^"M^/..a'/3' = 0. (C5) 

Here we have defined the differential operator 

= U8>i^8l + 2C^^ii\ (C6) 

and we have made use of the identity Q 

□ct = 4- (InA)'^cr;^, (C7) 



where A is the Van Vlcck-Morettc determinant (jB29ll . 

The solution to the differential equation (|G3|I for [/"/^" 1^ which is appropriate for the source term on the right 
hand side of Eq. (|G1|I is 

Ua0a'P' = A^/^ [ffQ'(Q5/3)/3' + KQa^ga'p'] , (C8) 



where we have used the identity l|B19|l . 

Next, we assume formal power series expansions for Vapa'p' and WafSa'p' of the form 



V^p^,p,ix,x') = V:f,^,p,{x,x')a". (C9) 



71=0 

and 



Wa.p^,p,ix,x') = J2 (CIO) 
n=0 

Note that these expansions do not define unique representations of the bitensors Vapa'p' and WafSa'p', since the 
coefficients VJ'^^/^/ and W^j^^,p, can be arbitrary functions of x and x'. However, one can obtain a unique set 
of coefficients from the following prescription Pick a bisolution W^^^ip, of the homogeneous wave equation 
'^ap^''W°^^,p, = 0. Then, substitute the expansions l|n9|l and HClOp into Eqs. ((^4)1 and (jUsl, simphfy using the 
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identities (jBSj) and ljC7p . and equate to zero the coefficients of powers of a. The result is the following recursive set 
of ordinary differential equations along the geodesic joining x and x' that allow one to solve for the coefficients: 



^afia'13' + 



n+1 



2n 



(Cll) 
(C12) 



and 



1 



2n(7i + 1 
^ -V" 



Equations (|C12|) and ljC13|) apply for n > 1. The power series (jC9|l and (|C10|I with these coefficients converge in a 
neighborhood of the diagonal x — x' for analytic metrics 1^. 

In this paper we are only interested in the coefficients V^p^,^,^ which can be obtained from Eqs. HC11() and (|C12|I . 
We expand each of these coefficients as covariant Taylor series of the form (|B12|) 



V-n 7 (5 

af3a'l3' — 9a' 9f3' 



(C14) 



where the coefficients v^. ^ are local tensors at x. We now specialize to the case n = 0, and substitute the expansion 
()C14(l for V^fj^ip, and the expression ljC8|l for Uapa'P' into the differential equation (|C11|I . We simpHfy using the 
identity (|B19p and the definition ljC6p . and expand the various terms as power series in cT;^ using the expansions 
(|B15|) . ljB25ll . (|B27|) . (|B32|) and IjBSSp . Equating the coefficients of the various powers of cT;^ then gives a series of 
equations that can be solved for the coefficients u° ^. The results are 



1)0 

"q/37(5 



71° 



"ap-fSeC 



Ca'jPSt 



2 ^afPS:€^ 



(C15) 
(C16) 

(C17) 



and 



V 



p 

apfSeCv 



— CajpS;6Cri + ^ (C'q7 ^eC'/3i5/jC ) ;»/ ^ g 5q7 (C'/j'^^eC^^i,^ ) 



(C18) 



where 



n, 



:9ajgp5 



:9aS9pj + Kgap9-tS- 



(C19) 



In Eqs. ljC15(l - (|C18|I . the right hand sides are understood to be symmetrized on the index par (a/3), on the index 
pair (7(5), and on as many of the index triplet {e^i]) as are present. When k = —1/2, the formulae IjClSI) - l|C17p 
agree with Eqs. (A20) - (A22) of AFO specialized to the vacuum case. 

We will also need the first two of the coefficients in the expansion of V^i^^i^,, which we obtain from Eq. ljC12|) with 
n = 1. The hard part of the computation is evaluating the source term on the right hand side of this equation. Using 
the definition (jC14p and the expansions (|B15|) . ljB25p and ljB27|l we obtain 



'-^^aPiS + -^^ap-fSe ^apjSe 



(C20) 
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Now inserting suitably symmetrized versions of the formulae (|C15ll - (|C18|I for the expansion coefRcients ^ and 
using the identity HB4() gives 



57 95 



12 



15 11 



a'' + 0(3^), 



(C21) 



where the right hand side is understood to be symmetrized on the index pairs (a/3) and (jS). Next, we expand both 
sides of Eq. (|C12p with n = 1 as a power series in a-p to 0{s), using the formulae (jClSp . (jC16p and (jC21|) to evaluate 
the right hand side, and the expansion (|C14|I with n = 1 on the left hand side. Equating the coefficients of the various 
powers of a.^p then gives a series of equations that can be solved for the coefRcients w„ ,^ . The results are 



7,1 



-^DCaj/3S + -^Ca^ p'^ C^pSv + -^^^Caj^^ CpSpu ~ '^2^9 ajC p^'^ ^ C S p^ p + y20 (Jtivp<y (J^^^^ (C22) 



and 



— □(Ca7/35;e) — -^^^afi-fSC pypa-eC^^'"^ + —gaigpsCpvpa-eC^^'"^ ~ —C^^^'^ CpSpiy;i 



80' 



720 



18 



12 



(■yPS-.n ^'-^erfy''^aCPs'^ 



r^^da-fCv^PS 



C -.v 



~1^9psCa'yriC-pC'^^'' e " QQ^PS^a C^r^epX ~ '^^9psCa-/r)ClP^^'^'^ e ~ '^^P&C^ CanepX 
+ Y^9psC^^'^'^Cjp„p-e, 



(C23) 



where again there is implicit symmetrization on the index pairs (a/3) and (7(5). When k = —1/2, Eq. (|C22|I agrees 
with Eq. (A23) of AFO specialized to the vacuum case^. 



APPENDIX D: EXPANSION COEFFICIENTS 



In this appendix we list the expressions for the coefficients Va...ri which appear in the expansion (|3.5() . These 
expressions are obtained by substituting Eqs. (|(yl5(l - IjC'lSp and l|('22(l - (|(]23p . specialized to k = 0, into Eqs. H3.6|l 
- The results are 



VaP'jSe 



' 2 ^a^yPS-jS I 



(Dl) 



V, 



12 1 



5: 



apf36 

I 

2 



a^ai^e'-'pyp'^ + T^9a'rCp'' ^C, 



g>-'a7 6^^7p ~t" '^ya-i^p e'-'Spr 
}_/^ pT 



3-5 7 

— HCa-^pS + -^Cj'p CjpSr + T^C^^^'^ CpSpr — -^ga-^g pS^ pr puC'^ 



(D2) 



^ Note that the expressions 1C151 - IClSt and 1C22I - 1C23I are all traceless on the index pair (75), aside from the terms involving the 
tensor H^^^^. This means that performing a trace reversal on the index pair (75) is equivalent to changing the value of k from to 
— 1/2, in agreement with the discussion after Eq. ICH above. 
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and 

Va/3(751ekp) 



12 



24 



12 

-5/37 



1 

30 



1^ T 



12 

A 

s -.p 



24 



12 



1 

20 



1 

90 



^ bar a\\^ z , 



e p 



90 ^ 



+5ep 
1 

~12 



A ;t 

/3<T 



1 

270' 



1 

90" 



A 

7 (5 ;(T 



+5. 



crp 



^5a75/35C'rAMi/C'^'^^'';j + — □(Ca7/35:e) 



X T , ^ /^X ; 



24 o . g- 



72 



12 



540' 



T^^dPt^^-reX-pCJ^^ - -^gptCsrXpCJ/'^ + ^5/37C5QrA:pC^ + l^gPjCt 



180 



540' 



360' 



rXfj, 



(D3) 



The right hand sides of Eqs. IjDll) and (|D2|) are understood to be symmetrized over the index pairs (a/3) and (7(5). In 
Eq. I|D3(I . as the notation on the left hand side indicates, we have only computed the piece of Vai3-ySeap which is totaUy 
symmetric on the indices {-jSap). This is because only this piece of Vap-ySea-p is needed for computing the expansion 
coefficient fa \ from Eq. (|3.14|) above. The right hand side of Eq. (|D3p is understood to be symmetrized on the index 
pair (a/3) and on the indices (jSap). 



APPENDIX E: SOME USEFUL DETERMINANT IDENTITIES 

One occasionally encounters determinants of rank two tensors in relativity. The normal expressions for evaluating 
determinants do not lend themselves naturally to expression using standard tensor notation. However, it is possible 
to write the determinant of any square matrix in terms of powers of the matrix and their traces, both of which are 
easily expressed in tensor notation. We list the three simplest here for a matrix A of various sizes 



Dimension 


Identity 










2x2 


DetA = 


[{TrA)^ - 


TriA^)] /2 




(El) 


3x3 


DetA = 


[{TrA)^- 


3{TrA)Tr{A^) H 


^2Tr{A^)] /6 


(E2) 


4x4 


DetA = 


[{TrAr- 


6{TrA)^ TriA^) 


+ 8{TrA) Tr{A^) + 3{Tr{A^))^ - 6Tr{A*)] /24 


(E3) 



Similar identities can be obtained for any dimension through a straightforward recursive application of Newton's 
Identities. 

To make the usefulness of these identities more apparent, we express the determinant in the 4x4 case in tensor 
notation using ljE3ll : 



DetA" 13 = 



(A"„)^ - 6 (A"„)' ^^7^-^/3 + SA'^^A'^^A'^sA^p + 3 (A'^pA^^)^ ~ eA'^^A^^A'^sA^J /24. (E4) 



While this expression is not a computationally efficient way of calculating the determinant, it clearly has the advantage 
of reducing the determinant operation to more familiar tensor operations. We include it here because we have not 
been successful in locating it in the literature. 
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